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PREFACE 

This  report  was  prepared  by  Darryl  J.  Calkins,  Research  Hydraulic  Engineer,  of  the  Ice 
Engineering  Research  Branch,  Experimental  Engineering  Division,  U.S.Army  Cold  Regions 
Research  and  Engineering  Laboratory. 

This  study  was  funded  under  DA  Program,  Ice  Engineering-,  Subprogram,  ice  Formation; 
CWIS  31350,  Format ional  Processes  of  Frazil  Ice. 

Dr.  D.  S.  Sodhi  and  Dr.  G.  D.  Ashton  of  CRREL  reviewed  the  technical  content  of  this 
report. 


ACCELERATED  ICE  GROWTH 
IN  RIVERS 

Darryl  J.  Calkins 

INTRODUCTION 

The  growth  of  an  ice  cover,  once  it  is  formed, 
results  from  the  continuing  heat  loss  from  the  cover 
to  the  atmosphere.  The  thickening  of  the  cover  is 
often  related  to  accumulated  degree-days  of  freezing 
with  an  appropriate  correction  factor  a for  the 
geographic  setting  or  heavy  snow  conditions.  The 
prediction  of  the  so-called  static  growth  of  the  ice 
cover  thickness  can  take  the  form  rjj  = a ( A 77 
where  r)i  = ice  thickness,  a = 0.6  to  0.9,  and  A 77  = 
the  accumulated  degree-days  below  freezing. 

The  effects  of  frazil  ice  deposits  beneath  an  ice 
cover  on  the  thickening  of  the  solid  ice  cover  have 
not  been  thoroughly  analyzed.  But  Bengtsson 
(1 978)  reports  ice  thicknesses,  due  to  frazil  ice,  up 
to  100%  greater  than  those  formed  by  normal  static 
growth  in  a river  in  Sweden.  This  report  considers 
the  impact  on  the  growth  of  the  solid  ice  sheet  due  to 
frazil  deposits  beneath  the  solid  ice  cover.  It  compares 
the  numerical  scheme  suggested  in  this  study  with 
field  data  collected  on  a river  having  approximately 
1 m of  frazil  beneath  a solid  ice  sheet. 

ANALYSIS  OF  HEAT  LOSS 

Figure  1 gives  a definition  sketch  of  the  problem. 

The  heat  flux  from  the  moving  water  to  the  frazil 
ice  deposition  is  given  as 

^wf  ~ MV^wf)  (I) 

where  7tw f is  a heat  transfer  coefficient,  Tw  is  the 
water  temperature,  and  Twf  is  the  melting  point 
of  the  ice  water  mixture,  0°  C.  Since  the  frazil 
ice  deposition  is  maintained  at  0°C,  the  temperature 
at  the  solid  ice/frazil  slush  interface  7"fj  is  also  0°C. 
Therefore,  there  is  no  heat  transfer  through  the 
frazil  slush  to  the  solid  ice;  and  until  the  frazil  ice 


is  removed  in  some  manner,  there  will  be  no  heat 
transfer  from  the  river  water  mass  to  the  solid  ice 
underside.  Consequently,  this  convection  heat 
term  can  be  neglected  in  the  analysis. 

The  thickening  of  the  ice  cover  is  then  related 
to  the  amount  of  heat  extracted  by  atmospheric 
sources.  The  heat  flux  0;  and  0S  through  ice  and 
snow  layers,  assuming  a linear  temperature  dis- 
tribution through  the  respective  thicknesses,  can 
be  represented,  respectively,  as 

*.  !'•  t 


0S  = *S  r's~r*a 

Vi 


where  /r(  and  are  the  thermal  conductivities  of 
the  ice  and  snow;  tjj  and  t]f  are  the  thicknesses 
of  the  ice  and  snow  layers;  Tfl  is  the  frazil  ice 
temperature;  Tjs  is  the  ice/snow  interface  tempera- 
ture; and  Tsa  is  the  snow  surface  temperature, 
which  is  < 0°C. 

The  heat  flux  from  the  snow  surface  (neglecting 
evaporation  or  condensation)  can  be  calculated  in 
a similar  manner  to  0wf  through  the  use  of  a heat 
transfer  coefficient  /?a$,  such  that 


h*  = hi%  (rsa-r») 


where  7"a  is  the  air  temperature.  Since  <t>,  = 0,  = 0S 
the  snow/ice  and  snow/air  interface  temperatures 
can  be  eliminated  in  eqs  2,  3,  and  4,  such  that 
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Figure  1.  Definition  sketch. 

ICE  COVER  GROWTH 

The  thickening  of  the  ice  cover  is  governed  by 
the  energy  balance  at  the  solid  ice/frazil  slush 
interface.  Since  the  heat  flux  to  the  solid  ice  from 
the  river  water  is  zero,  the  governing  equation  is 


(6) 


where  Pj  is  the  ice  density,  X is  the  latent  heat  of 
solidification,  and  dv-Jdt  is  the  rate  of  ice  thicken- 
ing. The  heat  of  solidification  is  the  energy  required 
to  change  the  phase  of  water  from  liquid  to  solid  at 
0°C. 

Frazil  ice  slush  comprises  a mixture  of  individual 
ice  crystals,  mostly  flocculated  frazil  discs,  and 
water.  The  growth  of  solid  ice  in  the  presence  of 
these  frazil  discs  can  increase  substantially  because 
the  frazil  ice  crystals  freeze  into  the  interstitial 
water  due  to  heat  loss  from  the  atmosphere.  Since 
the  frazil  ice  occupies  a given  percentage  of  a fixed 
volume,  the  amount  of  heat  necessary  to  convert 
the  ice-water  mixture  to  ice  is  correspondingly 
reduced  by  the  ratio  of  volume  of  water  to  the 
total  volume  of  the  mixture. 

A correction  for  the  latent  heat  of  solidification 
for  ice-water  mixtures  should  be  introduced.  The 
correction  is  related  to  the  porosity  <t>  of  the  frazil 
ice  suspended  in  the  fluid,  where  the  porosity  is 
defined  as  the  ratio  of  void  volume  (the  water 
fraction)  to  the  total  volume.  Consequently,  eq  6 
should  be  more  appropriately  written  as 


0.  - p,  0 X-~ 


and  combining  it  with  eq  5 yields 


(7) 


Tfrh 


(fljk,)  + (rjs/*s)  + 0//»a 


= Pi  <P  X 


dt 


(8) 


which  is  very  similar  to  the  equations  proposed  by 
Assur  and  Weeks  (1963)  and  Ashton  (1978). 

Because  of  the  dependence  of  some  parameters 
in  the  above  equation  on  time  (i.e.,  7"a,  »}s  and  /tas)( 
it  is  impractical  to  present  the  exact  solution  as 
one  would  be  continually  updating  these  parameters. 
An  analytical  integration  is  suggested  such  that 


Ar?i  = 


Tu-T, 


p,0X  (r},/k)  + (vjks)  + (1//»a 


At 


(9) 


NUMERICAL  RESULTS 

The  prediction  of  ice  thickness  can  be  written  as 

ijf*1  =uf+Ai?f  (10) 

where  the  superscript  k is  the  increment  of  time. 

The  values  of  various  parameters  have  been  assumed 
as  follows:  X = 3.34x  1 05  ) /kg,  p,  = 91 6 kg/m3, 

7sm  = 0°C,  k,  = 2.24  W/m-deg  and  At  = 1 day. 
Numerical  schemes  with  idealized  but  realistic  values 
for  given  meteorological  conditions  are  given  in 
Figures  2 and  3,  and  the  effect  of  frazil  ice  slush 
on  the  ice  covering  thickening  is  quite  apparent. 
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Figure  2.  Ice  growth  for  different 
frazil  ice  porosities. 
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Figure  3.  Ice  thickness  influenced  by  a 
combination  of  snow  and  frazil. 


Figure  4.  Ice  thickness  growth  at  Ottauquechee 
River,  Quechee,  Vermont. 


beneath  the  solid  cover.  Two  sections  were  chosen 
that  had  detailed  measurements  over  a 10-day  per- 
iod, with  minimal  snow  cover.  The  average  wind 
velocities  over  this  10-day  period  were  less  than 
4 m/s.  An  average  heat  transfer  coefficient  (25 
W/m2-°C)  suggested  by  Ashton  (1978)  for  open 
water  conditions  was  used;  and  an  average  thermal 
conductivity  of  the  snow  was  taken  as  0.3  W/m-deg 
after  Mellor  (1976).  The  ice  thickness  measure- 
ments were  resolved  to  0.0032  m.  The  snow  depth 
increment  was  taken  as  the  measured  value  at  a 
nearby  snow  course.  The  10-day  period  began 
with  no  snow  on  the  ice  cover  and  with  only  one 
storm  occurring  on  31  December,  with  a snowfall 
accumulation  of  roughly  0.05  m. 

The  numerical  scheme  was  run  for  two  dif- 
ferent sites  on  the  Ottauquechee  River,  Quechee, 
Vermont.  The  first  site  was  evaluated  using  0 = 1.0 
and  <t>  = 0.75;  at  the  end  of  10  days,  the  numeri- 
cal scheme  with  0 = 0.75  more  nearly  predicted  the 
final  growth,  although  the  predicted  value  was  low 
(Fig.  4).  When  the  ice  prediction  scheme  used 
<t>  = 1.0,  a significant  underprediction  of  roughly 
18%  was  noted.  For  the  second  site,  using  0 = 0.75, 
the  numerical  method  predicted  a slightly  greater 
ice  thickness.  The  insulating  effect  on  the  snow 
cover  is  very  evident  in  the  field  ice  thickness  data, 
but  the  numerical  method  does  not  reflect  the 
significant  decreased  growth  over  this  time  period. 
Since  neither  the  snow  thickness  nor  its  thermal 
conductivity  was  measured  at  the  time  on  the 
river  itself,  there  could  be  a substantial  difference 
in  values  due  to  just  snowdrifting. 

The  porosity  of  the  frazil  ice  was  not  measured 
but,  since  the  average  frazil  ice  depth  was  less  than 
1 m,  0 = 0.75  appears  a reasonable  value  at  this 
time. 


Dean  (1978)  made  field  measurements  on  the 
porosity  of  frazil  ice  slush  where  the  variations  of  DISCUSSION 

frazil  depth  were  from  2 to  1 2 m and  the  porosity 

values  were  0.48  to  0.34,  respectively.  Flowever,  The  prediction  of  ice  thickness  is  paramount 

more  information  from  the  field  is  needed  to  in  trying  to  establish  a safe  ice  sheet  to  support 

confirm  these  values  as  well  as  observations  of  heavy  loads.  The  growth  of  ice  cover  thickness 

any  changes  that  occur  throughout  the  winter  is  accelerated  by  frazil  ice  beneath  the  cover  due  to 

season.  the  reduction  in  the  effective  latent  heat  of  solidifi- 

cation. The  field  data  gathered  by  the  author  sup- 
port this  theory,  as  well  as  the  data  gathered  by 

FIELD  VERIFICATION  Bengtsson  (1978)  on  two  rivers  in  Sweden,  where 

excesses  of  0.1 2 m and  0.50  m were  measured  above 

To  verify  the  above  analysis,  field  data  were  the  so-called  static  growth  in  rivers  laden  with  frazil 

chosen  from  a small  river  with  frazil  ice  deposits  slush  beneath  the  solid  ice  cover. 


On  the  river  where  the  excess  ice  growth  ex- 
ceeded 0.50  m,  the  so-called  static  ice  growth  was 
roughly  0.50  m.  These  two  ice  thickness  measure- 
ments were  roughly  40  m downstream  of  each 
other  in  a rather  large  pool.  Just  upstream  of 
this  pool  was  a long  rapids  section  where  the 
frazil  ice  was  produced. 


CONCLUSIONS 

The  effects  of  frazil  ice  deposits  beneath  a 
solid  ice  cover  are  very  significant  in  increasing  the 
solid  ice  thickness.  The  effective  latent  heat  of 
ice  solidification  is  reduced  proportionally  to  the 
, porosity  of  the  slush  ice  beneath  the  cover.  The 

> increase  in  the  amount  of  solid  ice  thickness,  as 

shown  in  two  very  idealized  examples  in  this  study, 
can  be  roughly  50%  to  90%  depending  on  mete- 
orological conditions.  One  set  of  field  measure- 
ments from  Sweden  indicated  a 1 00%  increase  in 
solid  ice  thickness  caused  by  frazil  ice  slush. 

Ice  growth  rates  of  the  numerical  model  used 
and  field  data  were  compared;  and  the  numerical 
model  predicted  the  ice  thickness  after  10  days  to 
within  a reasonable  accuracy,  with  less  than  13% 
error  for  two  sites  studied.  Although  a more 
thorough  field  investigation  is  necessary  to  sub- 
stantiate some  of  the  physical  constants  of  snow 
and  ice,  the  results  of  this  study  clearly  indicate 
the  effect  of  frazil  ice  on  solid  ice  growth  rates. 
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